APPLICATION FOR UNITED STATES LETTERS PATENT 



FOR 



A FILTER FOR A MODULATOR AND METHODS THEREOF 



Inventors: Daniel YELLIN 
KobbyPICK 



Intel Reference No.: P17450 

EPLC Reference: P-6063-US 

Prepared By: Miriam Paton 

Eitan, Pearl, Lafzer & Cohen Zedek, LLP. 



Attorney Docket No.: P-6063-US 
Assignee: Intel Corporation 



A FILTER FOR A MODULATOR AND METHODS THEREOF 



BACKGROUND OF THE INVENTION 

[0001] In polar modulation, a signal is separated into its instantaneous amplitude and 
phase/frequency components (rather tiian into the classical in-phase (I) and quadrature 
(Q) components), and the amplitude component and phase/frequency component are 
modulated independently. The amplitude component may be modulated with any 
suitable amplitude modulation (AM) technique, while the phase/frequency component 
may be modulated using an analog phase locked loop (PLL). 

[0002] To allow reasonable operation, the bandwidth of the PLL may be quite small, 
much smaller than the actual bandwidth of the transmission signal's instantaneous 
phase/frequency. For example, in the case where the PLL is fed by a sigma-delta 
converter that has a high pass noise nature, the loop filter may be narrow enough to 
attenuate the sigma-delta quantization noise and the phase noise of the PLL. A pre- 
emphasis filter may emphasize, prior to modulation, those firequency components that 
would be attenuated by the PLL. The pre-emphasis filter may employ inverse 
filtering to the linearized response of tiie PLL. This inverse filtering may yield a 
high-order infinite impulse response (IIR), which may suffer from stability problems. 
[0003] Conventional practice involves calibration mechanisms in order to accurately 
calibrate the PLL to the predefined pre-emphasis filter. Without calibration, there is a 
risk that the pre-emphasis filter will not match the inverse to the PLL closed loop 
transfer fimction, which may result in enhancement of the phase noise. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0004] Embodiments of the invention are illustrated by way of example and not 

linodtation in the figures of the accompanying drawings, in which like reference 

numerals indicate corresponding, analogous or similar elements, and in which: 

[0005] FIG. 1 is a block-diagram illustration of an exemplary commimication system 

including a transmitter according to some embodiments of the invention; 

[0006] FIG. 2 is a block-diagram illustration of an exemplary fi:actional-N sigma- 

delta modulator according to some embodiments of the invention; 

[0007] FIG. 3 is a flowchart illustration of an exemplary method to determine a 

transfer function of a filter, according to some embodiments of the invention; 

[0008] FIG. 4 is a block-diagram illustration of an exemplary &actional-N sigma- 

delta modulator including an adaptive filter, according to some embodiments of the 

invention; 

[0009] FIG. S is a block-diagram illustration of an exemplary linearized 
approximation to the exemplary firactional-N sigma-delta modulator of FIG. 2, 
according to some embodiments of the invention; 

[0010] FIG. 6 is a block diagram of an exemplary system including an equalizer. 
[0011] It will be appreciated that for simplicity and clarity of illustration, elements 
shown in the figures have not necessarily been drawn to scale. For example, the 
dimensions of some of the elements may be exaggerated relative to other elements for 
clarity. 
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DETAILED DESCRIPTION OF THE INVENTION 

[0012] In the following detailed description, numerous specific details are set forth in 
order to provide a thorough understanding of embodiments of the invention. However 
it will be understood by those of ordinary skill in the art that embodiments of the 
invention may be practiced without these specific details. In other instances, well- 
known methods, procedures, components and circuits have not been described in 
detail so as not to obscure the description of embodiments of the invention. 
[0013] FIG. 1 is a block-diagram illustration of an exemplary communication system 
including a transmitter according to an embodiment of the invention. A 
communication device 100 is able to commimicate with a communication device 102 
over a commxmication channel 104. Although the scope of the invention is not 
limited in this respect, the communication system shown in FIG. 1 may be part of a 
cellular communication system (with one of communication devices 100, 102 being a 
base station and the other a mobile station, or with both communication devices 100, 
102 being mobile stations), a pager communication system, a personal digital assistant 
and a server, etc. A non-exhaustive list of examples for the communication system 
shown in FIG. 1 includes a Global System for Mobile Communications (GSM) 
system, a General Packet Radio Service (GPRS) system, an Enhanced Data for GSM 
Evolution (EDGE) system, code division multiple access (CDMA), CDMA2000, 
Wideband-CDMA (WCDMA), AMPS and any other wireless standard. 
[0014] Communication device 100 may include at least a transmitter 106 and an 
antenna 108. Communication device 102 may include at least a receiver 110 and an 
antenna 112. Antennas 108 and 112 may be of any desired kind such as, but not 
limited to, dipole, Yagi and multi-pole and the like. Moreover, communication device 
100 may include a receiver (not shown). Similarly, communication device 102 may 
include a transmitter (not shown), 

[0015] Transmitter 106 may include at least a baseband symbol generator 114 to 

generate a signal of baseband symbols, a splitter 116 to split the signal into its 

instantaneous amplitude and phase/firequency components, an amplitude path 118 to 

modulate and amplify the amplitude components, a fim:tional-N sigma-delta 

modulator 120 to modulate and up-convert the phase/frequency components, and a 

power amplifier 122 to amplify the output of fiactional-N sigma-delta modulator 120 
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Avith a gain controlled by the output of amplitude path 118. Antenna 108 is coupled to 
power amplijBer 122 to transmit the output of power amplifier 122. 
[0016] Baseband symbol generator 114 may be implemented in accordance with a 
wireless standard. Splitter 116 may be implemented in hardware, software or 
firmware or any combination thereof. 

[0017] FIG. 2 is a block-diagram illustration of an exemplary firactional-N sigma- 
delta modulator according to an embodiment of the invention, such as, for example, 
fiactional-N sigma-delta modulator 120 of FIG. 1. Fractional-N sigma-delta 
modulator 120 has a digital portion and an analog portion. The digital portion may 
include at least a differentiator 202, a filter 204 (such as, for example, a pre-emphasis 
filter), an amplifier 206, a summer 208, and a sigma-delta converter 210. The analog 
portion, which implements an analog phase locked loop (PLL) 212, may include at 
least a non-linear frequency divider 214, a reference oscillator 216 to produce a signal 
having a frequency , a non-linear phase-frequency detector 218, a loop filter 220, 

an amplifier 222 having a gain Ky^ and a voltage-controlled oscillator (VCO) 224. 

The input of VCO 224 is denoted Y. PLL 212 is also known as a "fractional-N phase 
locked loop" unit. Amplifier 222 and VCO 224 may be implemented as a single unit 
(VCO 223); the separation into two elements is for the sake of analysis. 
[0018] Reference oscillator 216 may produce a signal having a reference frequency 
Fj^j, . An example of reference frequency is approximately 26 MHz, although 
other reference frequencies may be used instead. 

[0019] Differentiator 202 may differentiate the phase symbols received from splitter 
116 of FIG. 1 to obtain the instantaneous frequency W of the baseband signal. 
Differentiator 202 may be implemented in hardware, software or firmware or any 
combination thereof. 

[0020] Filter 204 may be determined so that the instantaneous frequency W is 
transferred to the input 7 of VCO 224 with little or no distortion, as will be described 
in more detail hereinbelow. Any implementation of a digital filter is smtable for filter 
204. 

[0021] Amplifier 206 may normalize the output of filter 204, denoted X, to PLL 
reference frequency units. 
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[0022] Summer 208 may add the output of amplifier 206 with a number p . 
Both the integer number N and the non-integer number fi having a value between 0 
and 1 may be set according to an instruction from a base station regarding the average 
output frequency of a transmitter of a mobile station. For example, if reference 
frequency is approximately 26 MHz, then Nmdcy have a value in the range of 29 
- 32 so that N • F^^^ has a value of approximately 800 MHz. 

[0023] Sigma-delta converter 210 may convert the output of summer 208 into an 
integer number that represents the instantaneous frequency division ratio of the PLL. 
[0024] Non-linear frequency divider 214 may divide tiie output of VCO 224 by the 
integer number provided by sigma-delta converter 210 to produce a divided-frequency 
signal. 

[0025] Non-linear phase-frequency detector 218 may compare tiie divided-frequency 
signal to the reference frequency signal produced by reference oscillator 216. Non- 
linear phase-frequency detector 218 may produce a control signal that corresponds to 
the phase difference and/or frequency difference between the two signals. Any 
implementation of a phase-frequency detector is suitable for phase-frequency detector 
218. 

[0026] The control signal, after smoothing by loop filter 220 and amplification by 
amplifier 222, may be applied to VCO 224 so that VCO 224 synthesizes a modulated 
output carrier signal. 

[0027] As mentioned above, filter 204 may be determined so that the instantaneous 
frequency Wis transferred to the input 7 of VCO 224 witii litde or no distortion (e.g. 
the overall response from Wto Yis close to flat in the frequencies of interest). 
[0028] FIG. 3 is a flowchart illustration of an exemplary method to determine a 
transfer function C of filter 204, according to an embodiment of the invention. The 
method may include the following: 

a) building a linear model for PLL 212 and calculating the transfer ftmction H(w) 
from the ou^utA" of filter 204 to the input 7 to VCO 224 (block 302); 

b) optionally adding various impairments of the different PLL components (e.g. phase 
noises, variations of the PLL parameters from their nominal values, and the like) to 
the model (block 304) — the variations of PLL parameters from their nominal values. 
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such as, for example, capacitors, resistors, open loop gain, etc., may be due to 
production inaccuracies; 

c) deciding on a topology for the transfer function C of jBlter 204 (block 306), and 

d) calculating a transfer function C to minimize a predefined cost function related to 
the instantaneous firequency Wand the input 7to VCO 224 (block 308). 

[0029] In some embodiments, the transfer function C is a stable transfer function. 
[0030] Mathematically, the transfer function C may be expressed as follows: 
C = ArgMiric {Cost{W, Y)} . 

[0031] In one embodiment, the cost function may be the mean square error (MSE) 
between the instantaneous firequency Wand the iuput Tto VCO 224, as follows: 

Cost(^^Y)= E)^(t)—Y{t)^^, where f is a time variable and £ is an expectation 

operator. In yet another embodiment, the expectation operator is replaced by time 
averaging and the cost fimction becomes: 

^ r/2 

Cost(fV,Y) = ±^ ^W(t)-Y(tfdt 

where 7* is a design parameter determining the integration time window. 

[0032] In another embodiment, the cost function may be a weighted MSE in the 

frequency domain, as follows: 



Cost(W,Y)=^E 



jF{w}w{w)-Y(wfdw 



where P(w) is a user-defined, positive, weight function. For example, weight fimction 
P(w) may give more weight to those frequencies where spectral cleanliness is more 
important. 

[0033] In other embodiments, other cost functions may be used, such as, for example, 
a cost function that measures the spectral cleanliness of the overall transmission 
signal. 

[0034] In the event that an MSE or weighted MSE cost fimction is used, then block 
304 of the method may be redimdant if all impairments may be represented as 
additive noise terms, since different choices of transfer function C will not affect that 
total contribution of these additive impairments to the MSE or weighted MSE cost. 
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[0035] Certain cost functions may be minimized using equalization theory, as will be 
described hereinbelow for a particular example. 

[0036] Any suitable topology for transfer function C may be used in block 306. For 
example, transfer function C may be a jfinite impulse response (FIR) of order p (in 
which case stability is guaranteed), or an infinite impulse response (IIR) having a 
rational transfer function of orders (p^g)^ etc. 

[0037] FIG. 4 is a block diagram of an exemplary fractional-N sigma-delta modulator 
420 including an adaptive filter, according to some embodiments of the invention, 
Fractional-N sigma-delta modulator 420 is similar to firactional-N sigma-delta 
modulator 120 of FIG. 2, and therefore similar components are referenced with the 
same reference numerals and will not be described in further detail. Fractional-N 
sigma-delta modulator 420 may include at least an adaptive filter 404 (such as, for 
example, an adaptive pre-emphasis filter). Adaptive filter 404 may include at least 
filter 204, an analog-to-digital (A/D) converter 402 and an adaptive algorithm 403. 
[0038] Adaptive algorithm 403 compares the input to filter 204 (the instantaneous 
firequency W) to the input to VCO 223 (after digitization). Adaptive algorithm 403 
adapts filter 204 according to the comparison. For example, if the error between 
instantaneous firequency W and input Y is defined by a particular cost function, then 
adaptive algorithm 403 may reduce the error iteratively. A non-exhaustive list of 
examples of adaptive mechanisms include Least Mean Squares (LMS), Recursive 
Least Squares (RLS), and the like. Thus, any impairments, offsets, drifts, etc. of the 
analog portion of the PLL may drive the filter values to those values that minimize a 
pre-specified adaptive mechanism cost function such as, for example, a mean squared 
error cost function. Variations in the PLL, such as for example, variations in 
temperature, voltage, aging, etc., may be compensated for by the adaptive algorithm. 
[0039] Since adaptive algorithm 403 enables the digital values of filter 204 to be 
adapted to variations in the analog PLL, it may not be necessary to calibrate PLL 212 
to predefined values for filter 204. 

[0040] Block 302 of the method of FIG. 3 includes building a linear model of a PLL. 
FIG. 5 is a block diagram of an exemplary linearized approximation to the exemplary 
fi:actional-N sigma-delta modulator of FIG. 2. In this approximation, which is made 
in the phase domain, PLL 212 has been approximated by a linearized PLL 512, and 



7 



Attorney Docket No.: P-6063-US 
Assignee: Intel Corporation 

sigma-delta converter 210 has been approximated by a linearized sigma-delta 
converter 510 including an all pass filter and an additive noise w(Z - A) . 
[0041] Non-linear jfrequency divider 214 has been approximated by a linearized 
frequency divider 514 including an amplifier 511 having a gain of InFj^^y a 
diSerentiator block 530 having a transfer function of s, where s is a Laplace domain 
variable, a subtraction block 532, an amplifier 534 having a gain of l/(iV^ + and an 
integrator block 536 having a transfer fimction of 1/s. 

[0042] VCO 224 has been approximated by a block 524 having a transfer function of 
1/s. Non-linear phase-frequency detector 218 has been approximated by a subtraction 
block 518, and loop filter 220 has been approximated by a block 520 having a transfer 
fimction T(s). 

[0043] The transfer fimction from the outputs of filter 204 to input Tto block 524 is 
then given by the following expression: 

X{s) l + iKy/{N^/3))T{s)/s ^ 

[0044] The design of the transfer fimction C of filter 204 noiay then be performed as 
follows. In order that the instantaneous frequency Whe transferred to the input Y of 
VCO 224 with little or no frequency and phase distortion, the overall transfer fimction 

from Wto Y may be of substantially 0 dB gain up to a cut-off frequency fo and may 
have substantially linear phase up to the cut-off frequency fg. 

[0045] The cut-off frequency /q may be defined experimentally by observing the 
spectrum of the instantaneous frequency W signal. Alternatively, the cut-off 
frequency fg may be defined as narrow as possible while still meeting the 
requirements of a communication standard. Other definitions of the cut-off frequency 
fo are also within the scope of the invention. 

[0046] One straightforward solution is to define the transfer fimction C from WtoX 
up to the cut-off frequency fo as the inverse IIR to the transfer fimction Y(s)/X(s)^ as 
given by the following expression: 

W{F) " {K,l{N^P)\T{s)ls ' - - ^ • / < /o (Equation 2) 
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where j is the square root of -1 . 

[0047] However, if flie transfer function C were implemented as the inverse UR to the 
transfer function Y(s)/X(s)^ then the following problems might arise: 

a) zeros and poles may need to be added to the inverse IIR to stabilize the filter; 

b) additional poles or filters may need to be added to the inverse IIR to attenuate the 

firequencies above the cut-off firequency fo so that the sigma-delta converter will not 

be forced into saturation and so as not to increase the quantization noise; and 

c) the order of the transfer function C would no longer be a design parameter; radier, 
it would have a one-to-one relation to the order of the closed loop transfer function. 
Therefore, according to some embodiments of the present invention, the transfer 
function C is not the inverse IIR to the transfer function Y(s)/X(s). Rather, the transfer 
function C of the filter is calculated to optimize predefined optimization criteria, as 
described hereinabove with respect to FIG. 3. 

[0048] Block 302 of the method of FIG. 3 includes calculating the transfer function C 
to minimize a predefined cost function related to the instantaneous firequency W and 
the input Fto the VCO of the PLL. The transfer function C may be calculated in a 
manner similar to the calculation of a minimum mean squared error (MMSE) 
equalizer, as will be explained in greater detail. FIG. 6 is a block diagram of an 
exemplary system including an equalizer. 

[0049] A known digital input signal u(i) may be subjected to an impulse response /t, 
where die impulse response h is defined as the bi-linear transform of the transfer 
function Y(s)/X(s). 

[0050] An additive white Gaussian noise n^WGN(i) may be shaped by a shaping 
filter g. Noise shaping filter g may be a high pass filter with a cut-off above 

frequency /q^ so that the equalizer response at high frequencies will be attenuated. 

[0051] The combination of shaped noise n(i) and the output of impulse response h is 
denoted v(i). The digital signal v(i) may be subjected to a FIR filter C to yield a 
digital signal i}(z). 

[0052] FIR filter C may be calculated analytically to minimize the mean square error 
that is defined as 

e = - . (Eqxaation 3) 
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[0053] The followiag e^qsressions may be useful in the calculation of FIR filter C: 



V=HU+N, 
'h{M + L) 



h{L) ••• h{-M + L) 



Us 



N = 



h{M) h(0) 

h{M-L) 
u{i-M) 

u{i) 
»(z + Z)' 

«0-4 



hirM) 



h{-L) h{-M-L) 



(Equation 4a) 
(Equation 4b) 

(Equation 4c) 



[(2£+lH2W+0] 



(Equation 4d) 



(Equation 4e) 



J[(2i+l>l] 

-vv^ore the bar denotes vector notation, the superscript H denotes the conjugate 
transpose, 2M+1 is assumed to be the length (i.e. number of coe£&cients) of impulse 
response h, and 2L+1 is assumed to be the length (i.e. number of coefficients) of the 
equalizer. 

[0054] The following assumptions may be made: 

1. There is no correlation between the error and .the observations. Thus, 

}= 0^ , "w^ere, again, £ is an expectation operator. 

2. There is no correlation between the input signal and the shaped noise. Thus 

= 0, VA: 5 where the asterisk denotes the complex conjugate. 

3. The input signal is independently distributed (ID). Thus 
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4. n AWGN is additive white Gaussian noise. 



,Jfc = 0 



Thus 
where 



E^_AWGN{iyn_AWGN{i + k)'}= ^^^^ 

<jI_awgn = ^{ \n_AWGN{}^ } is a fixed value for all i. 
[OOSS] From assumption 1 it follows that 

E^{})-V"}= E^i^yV"). (Equation 5) 

[0056] By substituting Equations 4a, 4b, 4d and 4e into Equation 5, it follows that: 
E^{})-V"]=C" E^ V"] 

^C" -ElH-U + Nyp" -H" 



= C" 



e{h -U-U" H"]+ e{n • N"}]' 



(Equation 6) 



where for simplicity, o-y (f) is denoted by cr^ , and where matrix G is dejSned in terms 

of shaping filter in a similar manner to the definition of matrix H in terms of 
impulse response h. 

[0057] Similarly, using Equation 4b and assumptions 2 and 3, it follows that: 

7"}= E\u{i) -{h-u+n)} 

= <jl'[h{L) - h(p) ... h(rL)Y 
[0058] From Equations 5, 6 and 7, and applying the complex conjugate operator, the 
MMSE analytical calculation of the FIR filter C is as follows: 



(Equation 7) 



C =[h h" -tJu +G G" ■al .ijff^Y 



h{L) 
h{0) 



(Equation 8) 



[0059] FIR filter C may be calculated empirically. For example, a simulation 
environment may be built in wduch u(i) is the instantaneous firequency at flie 
baseband, h(i) is the total response firom the input to the sigma-delta converter until 
the input to the VCO, shaping filter g(i) is chosen so that shaped noise n(i) will be a 
high pass noise a cutoff firequency of which is beyond the band of interest of the total 
response. A simulation may then be run to detemiine empirical values for v(i). The 
following quantities may then be calculated: 
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(Equation 9a) 



(Equation 9b) 



[0060] The empirically calculated FIR MMSE filter is then given by: 




(Equation 9c) 



[0061] While certain features of the invention have been illustrated and described 
herein, many modifications, substitutions, changes, and equivalents will now occur to 
those of ordinary skill in the art. It is, therefore, to be understood that the appended 
claims are intended to cover all such modifications and changes as fall within the 
spirit of the invention. 
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